Explosive condensation in a mass transport model 
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We study a far-from-equilibrium system of interacting particles, hopping between sites of a Id 
lattice with a rate which increases with the number of particles at interacting sites. We find that 
clusters of particles, which initially spontaneously form in the system, begin to move at increasing 
speed as they gain particles. Ultimately, they produce a moving condensate which comprises a finite 
fraction of the mass in the system. We show that, in contrast with previously studied models of 
condensation, the relaxation time to steady state decreases as an inverse power of In L with system 
size L and that condensation is instantenous for L — > oo. 
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Recent studies in non-equilibrium statistical physics 
show that diverse phenomena such as jamming in traffic 
flow [2], polydisperse hard spheres wealth condensa- 
tion in macroeconomies Q, hub formation in complex 
networks 0| , pathological phases in quantum gravity [j| , 
and general problems of phase separation @ can be un- 
derstood by the condensation transition. Condensation 
occurs when the global density of a conserved quantity 
(mass, wealth etc.) exceeds a critical value, and manifests 
itself as a finite fraction of the total system mass localized 
in space. A well-studied, fundamental model is the Zero- 
range Process (ZRP) which may serve as either a micro- 
scopic or effective description of non-equilibrium conden- 
sation In this model particles hop to the right on 
a closed chain of L sites with rates u(m) = 1 + 7/771 de- 
pending only on the number of particles to > at the de- 
parture site. The condensate, which exists in this model 
for 7 > 2 when density of particles is above some criti- 
cal value, remains static once it has formed, melting and 
reforming very rarely Q . This is caused by attractive in- 
teractions between particles expressed in u(m): the more 
particles are in the condensate, the slower it evolves. 

In this work we demonstrate a novel mechanism of non- 
equilibrium condensation motivated by processes such 
as gravitational clustering [Io| , formation of droplets in 
clouds or on inclined surfaces due to collisions [111, and 
differential sedimentation (l2| , where aggregation of par- 
ticles speeds up in time as a result of increasing exchange 
rate of particles between growing clusters. For example, 
raindrops falling through the mist increase their veloc- 
ity when gaining mass, which causes them to accrete 
mass even faster. To better understand the difference 
between the dynamical nature of the condensate in such 
processes and the static condensation which has previ- 
ously been studied [2-0, [HI, we consider a microscopic 
model of particles hopping between sites of a Id lattice 
with rate u(m, n) ~ (mn) 7 which increases with the num- 
bers to, n of particles at interacting sites. We shall show 
that for 7 > 2 condensation occurs through a contrast- 



u(m, n) ~ (mn) 1 



Figure 1: Model definition: a particle hops from site with m 
to site with n particles with rate u(m,n). 



ing dynamical mechanism to that previously considered 
- the formation of the condensate happens on a very 
fast time scale and we term it explosive. By considering 
the microscopic processes of the dynamics we show that 
the condensate moves with speed v ~ L 1 which increases 
with system size L. We argue that each cluster of parti- 
cles has a chance to develop into the condensate in finite 
time. It then follows from extreme value statistics that 
the time to form of the condensate decreases with sys- 
tem size as (InL) 1-7 , in contrast to ZRP. This counter- 
intuitive result means that condensation is instantenous 
for L — > 00. 

Model definition: The model we consider comprises M 
particles hopping to the right between sites of a peri- 
odic chain of length L as in the ZRP. Although partial 
asymmetry may also be considered, we restrict ourselves 
here to the case of totally asymmetric hopping: a parti- 
cle hops from site i to site i + 1 with rate u(mi,m,i + i) 
where rm, TOi+i are the occupancies of the departure and 
arrival sites, respectively. We assume the factorized form 



u(m,n) = (v(m) — v(0))v(n), 



(1) 



where the function 7; (to) grows as a power of to 

7; (to) = (e + to) 7 ~ to 7 , (2) 

with e <C 1 and 7 > [19]. Equation (fTJ) implies that 
u(0,n) = and that for large m, n, u(m,n) ~ (mn) 1 is 
the bigger the more particles are located on both sites. 
This has dramatic consequences for the dynamics. Com- 
paring simulation results of this model to ZRP dynamics 
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Figure 2: Top: state of the system {m;} for v(m) = (m+0.1) J 
and L = 100, M — 400 at different times: (a) initial con- 
dition with randomly distributed particles, (b) the rise of 
microscopic clusters of particles, separated by empty sites, 
(c) the macroscopic cluster (condensate) close to the steady 
state. Middle: positions of five most occupied sites (red 
squares for the largest cluster) as a function of time in this 
model. Bottom: the same plot for ZRP condensation for 
u(m,n) = 1 + 3/m shows completely different dynamics. 



in Fig. [5] reveals some striking differences (see also anima- 
tions in Supp. Material 14]). In ZRP, initial microscopic 
clusters are first formed, but they coalesce and grow 
quickly, until two macroscopic clusters are left. These 
slowly merge into the final macroscopic condensate by ex- 
changing particles through the other sites which form the 
fluid background. In our model, particles also aggregate 
into clusters (see Fig. [5Jd) but then these clusters start to 
move in the direction of hopping particles. This process 
speeds up in time; some clusters move faster as they gain 
particles in collisions, and one of them - the condensate - 
starts to dominate (Fig. [2t). Due to the rapid nature of 
this process we call it explosive condensation. The speed 
v = di max /dt at which the condensate travels through 
the system stabilizes after the system reaches the steady 
state. The motion of the condensate is similar to the 
"slinky"- like motion of a non-Markovian model [ll|. Fi- 
nally, smaller clusters move in the opposite direction to 
the main condensate at each collision. 



The dynamics thus differs significantly from the zero- 
range process. Surprisingly, both models share simi- 



lar static properties. In fact, they belong to a class 
of processes that have the important property that the 
steady state probability P({rrii}) of a configuration with 
mi, . . . , ttil particles at sites i, . . . ,L factorizes: 



P({mi})=l[f(r 



i=l 



with f(n) defined as 



(3) 



(4) 



Equation ([3]) requires two conditions on u (m,n) 01, 
which are satisfied for our model ([T]) and the ZRP (where 
u(m, n) = 1 + j/m for to > and w(0, n) = 0). In both 
cases we can choose /(l) = /(0) and obtain from Eq. ((3J) 
the large m behaviour f(m) ~ m -7 . It is known Q that 
for a power-law f(rn), condensation happens when the 
density of particles p = M /L exceeds the critical density 
p c = linwzF'(z)/F(z), where F(z) = E,JH^ 
and z plays the role of fugacity. For 7 > 2, p c is finite but 
for 7 < 2, p c — > 00. Therefore, condensation is possible 
only for 7 > 2 and for p > p c , which marks the transition 
between condensation/no condensation regimes Q- 

We now come back to the dynamics of our process 
and investigate what determines the speed of clusters 
and the condensate, how the clusters collide, how long 
it takes to reach the steady state, and how this time 
depends on the initial condition. We are interested in 
the limit of large M,L and fixed density p = M/L. 
For our choice v(m) = (e + to) 7 and 7 > 2, we obtain 
/(to) = /(0)e 7 m -7 for to > and the critical density, 

Pc = Em m /H/Em/H ~ e7 C(7 where 
C(7 — 1) is the Riemann zeta function. As the critical 
density is low one can make the simplifying approxima- 
tion that the clusters move in an otherwise empty system. 

Let us first calculate the speed at which the cluster of 
to particles moves through the system. We assume that 
at t = the cluster occupies site i, so that to^ — to, and 
that there are no particles at sites i — l,i + 1. The time r 
it takes to move the cluster to site i + 1 is the sum of times 
t m , t m -\, . . . , t\ it takes to move one particle to the right 
if the cluster has m, m — 1, . . . , 1 particles, respectively. 
Each nth hop is a random process with average duration 
t n given by the inverse of the hopping rate u(n, m — n), 
thus 



^— ' u(n, m — n) 

n—l x ' 



(5) 



Recalling that for condensation we are interested in 7 > 2 
and using Eqs. (Q][2]), we obtain that (r) « (em) -7 , which 
shows that larger clusters move faster. The condensate, 
which has M particles, moves 1/ (r) ~ (eAf) 7 = 
(ep) 7 L 7 sites per unit time, in agreement with simu- 
lations: for parameters from Fig. [5] we have measured 
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Figure 3: Top (a): A collision of two condensates having ini- 
tially m,i and mi+2 particles (left) proceeds through exchange 
of particles at site i + 1 (middle pictures). Arrows of dif- 
ferent sizes show relative magnitudes of hopping rates. Af- 
ter the collision (right), the masses are m'i and m' i+2 , with 
(m'i — mi) < 0, (m' i+2 — m i+ 2) > 0. Bottom left (b): ex- 
ample of stochastic evolution of mi (t) (red), rrii+\(t) (black), 
and mi+2 (t) (blue). Bottom right (c): probability distribution 
P(Am) of the difference Am = rrii — for = 25, 50, 100 
(circles, squares, diamonds) and m,+2 = 10, e = 0.1,7 = 3. 
In all cases (Am) « 0.4. 



the speed 63850 ± 100 whereas the formula (eM) 7 gives 
64000. 

To understand what happens when two condensates 
collide with each other, we assume that a bigger con- 
densate with TO, particles approaches a smaller one with 
m i+2 particles from the left, and that they are separated 
by an empty site i + 1, see Fig. (3Ji. Initially, the dynam- 
ics is dominated by hops from site i to site i + 1 because 
u(rrii, m l+ i) oc m] is bigger than u(m i+ i,m 1+2 ) oc m] +2 . 
As particles accumulate at site i+1, u(rrii+i, mi+2) grows 
and u(rrii,mi+i) decreases until they become compara- 
ble. This happens when rm ~ mi+2 since u(m, n) ~ 
(mn) 7 is symmetric for large m, n. Then the second 
half of the process becomes a time-reversed and space- 
inverted version of itself (see Fig. \3]p) ; one can think of 
the flow from i to i + 2 as a flow from i + 2 to i in 
reversed time which is identical to the flow from i to 
% + 2 in the first half of the process. Due to this time- 
reversal symmetry, the final configuration becomes the 
initial configuration reflected around site i + 1, modulo 
random fluctuations. In the case of ((2]), the slightly bro- 
ken symmetry of u(m, n) produces a small net current of 
particles from smaller to bigger clusters (see Fig. [2k). 

We now venture to draw the following picture of con- 
densation dynamics. First, small clusters are formed ran- 
domly from the initial state. For a system of size L there 
will be TV = O(L) such clusters. Subsequently, these 
clusters move ballistically between collisions, which are 
almost elastic. One of them soon collects more parti- 
cles than the rest and starts moving at increasing speed, 
gaining mass and becoming the final condensate. Let 




Figure 4: Sketches of f(T) (left) and h(T ss ) = 
cxp(-iV f^"" f(T)dT) (right). h(T ss ) can be approximated 
by a step function (see text). 



us calculate the time T ss for the system to relax to sta- 
tionary state. Each cluster will go through a series of 
collisions and either dissolve into the background or be- 
come the condensate; in either case we can associate a 
relaxation time T to each cluster (with T — 00 if the 
cluster disappears). Then T ss will be the minimal time 
out of Ti, . . . , Tjv relaxation times for all TV clusters: 



T ss = mm{Ti, . . . ,T N } 



(6) 



The relaxation process of a particular cluster of initial 
mass mo is a series of transitions at times t n at which it 
moves by one site to the right and (possibly) exchanges 
a chunk of mass Am n with other clusters: 



t n -l ^ 



+- Am r , 
At n . 



(7) 

(8) 



Here At n is the time between two jumps and is exponen- 
tially distributed as 



p n (At n ) = \ n e 



-A„At„ 



(9) 



where A„ oc m 7 is the speed of the cluster. Let us cal- 
culate the probability distribution f(T) of the relaxation 
time T = At\ + At2 + . . .. Numerical simulations suggest 
that the mass m n increases linearly through the colli- 
sions. We may thus assume that m n oc n and A„ ~ An 1 
for large n. Then T is a sum of independent exponential 
random variables and f(T) is given by 



f(T) 



1 

2^ 



(10) 



where /(w) is the product of characteristic functions of 
exponential distributions ((9]): 



/>)= um»>)= n i 



1L 1 - iuj/\ n 



(11) 



We expect that f(T) has the shape depicted in Fig.@]and 
that it decays to zero for T -+ 0. The large-w behaviour 
of /(w), which corresponds to small-T behaviour of f(T), 
is given by [3| 



/M--z(27r) 7 / 2 W^exp 



2iri(iuj/A) 1 ^ 
e 2 W7 — l 



(12) 



Now. we must invert the Fourier transform to recover 

f(T). For small T, this may be done by the saddle point 

approximation to the integral over uj (dominated by uj = 
0(^-7/(7-1))) and Qne obtaing 



f(T) oc CT^-WMi-^ exp -B(AT) 



(13) 



where B, C are some real, positive constants. If we as- 
sume that each cluster evolves independently, the relax- 
ation time ^ of the system becomes the minimum out 
of N independent random variables distributed accord- 
ing to f(T). Extreme values statistics tells us that the 
distribution P(T SS ) is given by 



P(T SS ) = Nf(T ss ) 



f(T)dT 



N-l 



(14) 



and integrating by parts and expanding for f(T) small, 



(Tss) 



exp -N / f(T)dT dT si 



(15) 



Knowing the small- T behaviour (fT3")) of f(T), we can cal- 
culate the average (|15[) for large N as follows. The func- 
tion h(T ss ) = cxjp(—N J Q Tss f(T)dT) approaches a step 
function for large N , see Fig. [4j The integral (fl~5|) over 
T ss , then becomes (T ss ) = / °° h(T ss )dT ss = to, where t 
is the position of the step in h, which can be identified 
as the point at which h" = 0, yielding 



f'(t ) = Nf 2 (t ) 



(16) 



Inserting the short-time behaviour (]13fl of / into this con- 
dition one obtains 



CN(j - 1) = Bexp B(At 



A 7-1 i n 



Taking logarithms yields 



to 



In N - In 



^U-ilnio 

7-iy 2 



1-7 



(17) 



(18) 



Thus, recalling N = O(L) and (T ss ) = in, the relaxation 
time asymptotically decreases as 



{Tss) = c 2 (c 3 + InL) 1-7 . 



(19) 



This form crosses over from (T SB ) = l/(const + O(lnL)) 
for small L to (T ss ) = O^mi) 1 " 7 ) for large L. This dif- 
fers much from ZRP-like models where (T ss ) ~ L 2 grows 
with L 0, . Since the time to steady state decreases 
with L, an infinite system relaxes instantaneously. This 
is reminiscent of instantaneous gela tion known from the 
theory of coagulation processes [171 . In fact, our model 
provides a non-trivial example of instantaneous gelation 
in a spatially-extended system. However, the model and 
its effective description in terms of colliding clusters dif- 
fer from coagulation processes in that there is exchange 
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Figure 5: {T ss ) _1 obtained in numerical simulations (points) 
and from formula (|19[) fitted to data points (lines). In all 
cases the density p = 2 and 7 = 3, 4, 5 (curves from bottom 
to top). Left: v(m) = (0.3 + m) 7 , every 5th site has initially 
10 particles. Right: v(m) = (1 + m) 7 particles are distributed 
randomly in the initial state. (Tss)" 1 for different 7 differ by 
orders of magnitude and hence they have been rescaled to plot 
them together. 



of particles between clusters rather than coagulation (a 
model with exchange of particles has been studied in 
Ref. [l8j|, see Supp. Material for more details). 

In the above derivation we assumed that A„ is strictly 
proportional to n 7 , and that the proportionality coeffi- 
cient is the same for all clusters. This is valid only if all 
clusters have the same initial size m = 1. To account 
for fluctuations of cluster sizes one should take the prod- 
uct (|11[) not from n = 1 but from some uq > 0, with 
no changing from cluster to cluster. However, this does 
not modify the asymptotic behaviour of f(uj), it only 
increases the constant C3 in Eq. p^|) . We have checked 
numerically evaluating Eqs. (fill . (fTU|) and (fl"5j) that (T ss ) 
for no > 1 has much stronger finite-size corrections and 
behaves as ~ l/(const + O(lnL)) for a wide range of L. 
Although C3 may in principle be calculated from our the- 
ory for n > 1, in practice it is simplest to treat C3 as a 
free parameter. In this way Eq. (|19[) fits numerical simu- 
lations very well. To check this, we measured the time it 
took the biggest cluster to reach the mean steady-state 
size of the condensate, M — Lp c . In Fig. Owe compare 
Eq. (|19| with (T ss ) obtained in simulations, for different 
initial conditions. We plot (T ss ) 1 because it shows con- 
vincingly that (T ss ) 1 grows to infinity for L — > 00, and 
therefore (T ss ) — > in this limit. 

In conclusion, we have elucidated a form of dynamic 
condensation which happens in far-from-equilibrium sys- 
tem of hopping particles. In contrast to previously stud- 
ied models, the condensate moves through the system 
and its dynamics speeds up in time — hence we term the 
condensation "explosive". The relaxation is dominated 
by the process of initial coalescence which is the slowest 
stage of condensate formation, at variance with previ- 
ously studied models of condensation such as the ZRP 
where this stage is the fastest. It remains to be seen 
whether condensation can be made "explosive" also in 
models which do not have a factorized steady state, such 
as those with spatially extended condensates fljjj . 
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